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We simulate the central reactions of 20 Ne+ 20 Ne, 40 Ar+ 45 Sc, 58 Ni+ 58 Ni, 86 Kr+ 93 Nb, 
129 Xe+ 118 Sn, 86 Kr+ 197 Au, and 197 Au+ 197 Au at different incident energies for different equations 
of state, different binary cross sections and different widths of Gaussians. A rise-and-fall behavior 
of the multiplicity of intermediate mass fragments (IMFs) is observed. The system size dependence 
of peak center-of-mass energy E™ 3 ; and peak IMF multiplicity (Ni MF) max is also studied, where 
it is observed that E™^, x follows a linear behavior and (NiMF) max shows a power-law dependence. 
A comparison between two clusterization methods, the minimum spanning tree and the minimum 
spanning tree method with binding energy check (MSTB), is also made. We find that the MSTB 
method reduces the (N/mf)" 1 " 1 , especially in heavy systems. The power-law dependence is also 
observed for fragments of different sizes at E™^, 3 ; and the power-law parameter r is found to be close 
to unity in all cases except A ma:r . 
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I. INTRODUCTION 

At high excitation energies, the colliding nuclei may 
break into several small and intermediate size fragments 
and a large number of nucleons are also emitted [H-dj]. 
The emission of intermediate mass fragments (IMFs) in 
nuclear collisions was studied for more than a decade dur- 
ing which several experimental groups have carried out a 
complete study of fragment formation with 4ir detectors 
[fJ-UJj] . These studies revealed that the fragments formed 
in heavy-ion collisions depend crucially on the bombard- 
ing energy and impact parameter of the reaction 
Therefore, these experimental studies of fragmentation 
offer a unique opportunity to explore the mechanism be- 
hind the formation of the fragments. Moreover, one can 
also pin down the role of dynamics in fragment formation 
and its time scale. 

Recently, there has been increasing interest in the ef- 
fects of reaction dynamics on the production of IMFs 
and light charged particles (LCPs, Z = 1 or 2). Sisan 
et al. 6| studied the emission of IMFs from central col- 
lisions of nearly symmetric systems using a 47r-array set 
up, where they found that the multiplicity of IMFs shows 
a rise and fall with increase in the beam energy. They 
observed that E™^f (the energy at which the maximum 
production of IMFs occurs) increases linearly with the 
system mass, whereas a power-law (oc A T ) dependence 
was reported for peak multiplicity of IMFs with power 
factor t = 0.7. Peaslee et al. Q, however, studied the 
asymmetric system 84 Kr+ 197 Au in the incident energy 
range from 35 to 400 MeV/nucleon and obtained an en- 
ergy dependence of multifragmcntation. Their findings 



revealed that fragment production increases up to 100 
MeV/nucleon and then decreases with increase in inci- 
dent energy. De Souza et al. t 8ij studied the central col- 
lisions of 36 Ar+ 197 Au from 35 to 120 MeV/nucleon and 
observed that IMF multiplicity shows a steady increase 
with increase in the incident energy. The IMF multiplic- 
ity decreases, however, when one moves from central to 
peripheral collisions. However, Tsang et al. in their 
investigation of 197 Au+ 197 Au collisions at E/A = 100, 
250, and 400 MeV, found the occurrence of peak multi- 
plicity at lower energies for central collisions whereas it is 
shifted to higher energies for peripheral collisions. Stone 
et al. Q used a nearly symmetric system of 86 Kr+ 93 Nb 
from 35 to 95 MeV/nucleon to obtain IMF multiplic- 
ity distribution as a function of beam energy by select- 
ing central events. Ogilvie et al. [To[ also studied the 
multifragment decays of Au projectiles after collisions 
with C, Al, and Cu targets at the bombarding energy 
of 600 MeV/nucleon using the ALADIN forward spec- 
trometer at GSI, Darmstadt, with the beam accelerated 
by Schwerionensynchrotron (SIS). They found that, with 
increasing violence of collision, the mean multiplicity of 
IMFs originating from the projectile first increases to a 
maximum and then decreases again. As mentioned ear- 
lier, Sisan et al. 6] reported that the peak multiplicity 
of IMFs as well as peak center-of-mass energy scale with 
the size of the system. In a recent communication, Ver- 
mani and Puri [ll| succeeded partially in explaining the 
above-mentioned behavior by using the quantum molec- 
ular dynamics (QMD) approach. Here we plan to extend 
the above study by incorporating various model ingredi- 
ents such as equation of state, nucleon-nucleon (nn) cross 
section, and Gaussian width. The role of different clus- 
terization algorithms is also explored. We attempt to 
find out whether these ingredients have sizable effects. 
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II. THE FORMALISM 
A. Quantum molecular dynamics model 

We describe the time evolution of a heavy-ion reaction 
within the framework of a QMD model |l|-|3j, |l2[ , which 
is based on a molecular dynamics picture. Here each 
nucleon is represented by a coherent state of the form 



4> a (xi,t) 
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= _ p-(*l-*a(t))< 



gip a (a:i-a;a) e -- 



(1) 



Thus, the wave function has two time-dependent param- 
eters, x a and p a . The total n-body wave function is 
assumed to be a direct product of coherent states: 



<j) = 4> a (xi,X ot) p ot ,t)(pp(x2,Xp,P/3,t). 



(2) 



where antisymmetrization is neglected. One should, how- 
ever, keep in the mind that the Pauli principle, which is 
very important at low incident energies, was taken into 
account. The initial values of the parameters are cho- 
sen in a way that the ensemble (Ar+Ap) nucleons give 
a proper density distribution as well as a proper mo- 
mentum distribution of the projectile and target nuclei. 
The time evolution of the system is calculated using the 
generalized variational principle. We start out from the 
action 



S = C[4>, 0*1 ^t, 



with the Lagrange functional 



C = 

















dt 





(3) 



(4) 



where the total time derivative includes the derivatives 
with respect to the parameters. The time evolution is 
obtained by the requirement that the action is stationary 
under the allowed variation of the wave function 



5S = 5 



C[<j>, (j>*]dt = 0. 



(5) 
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If the true solution of the Schrodinger equation 
is contained in the restricted set of wave function 
ef> a (xi, x a ,p a ) , this variation of the action gives the ex- 
act solution of the Schrodinger equation. If the parame- 
ter space is too restricted, we obtain that wave function 
in the restricted parameter space, which comes close to 
the solution of the Schrodinger equation. Performing the 
variation with the test wave function (2), we obtain for 
each parameter A an Euler-Lagrange equation: 



d_dC dC_ {) 



(6) 



For each coherent state and a Hamiltonian of the form, 



the Lagrangian and the 



Euler-Lagrange function can be easily calculated [12[ : 
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Thus, the variational approach has reduced the n-body 
Schrodinger equation to a set of 6n different equations 
for the parameters that can be solved numerically. If 
one inspects the formalism carefully, one finds that the 
interaction potential, which is actually the Bruckner G- 
matrix, can be divided into two parts: (i) a real part 
and (ii) an imaginary part. The real part of the poten- 
tial acts like a potential, whereas the imaginary part is 
proportional to the cross section. 

In the present model, the interaction potential com- 
prises of the following terms: 



V a p = V? oc + V^ oc + Vcoui + V Y uk + V: 



MDI, 



(10) 



where Vi oc is the Skyrme force and Vcoui, Vy u k and Vmdi 
define, respectively, the Coulomb, Yukawa and momen- 
tum dependent potentials. The Yukawa term separates 
the surface which also plays a role in low- en ergy pro- 
cesses like fusion and cluster radioactivity [l3|, [l4J . The 
expectation value of these potentials is calculated as 



• I n 



v Inr 



/ u Pa ) , 

xd 3 r a d 3 rpd 3 pad 3 pp, (11) 
i)//3(P/3,ra,i)/ 7 (p 7 ,r 7 ,i) 



xV I {3) (r a ,rp,r J )d 3 r a d 3 rpd 3 r 1 
xd 3 p Q d 3 p ( gd 3 p 7 . 
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where f a (p a ,r a ,i) is the Wigner density which corre- 
sponds to the wave functions [Eq. 2] . If we deal with the 
local Skyrme force only, we get 



V 



Skyrme 



A T +A f 

= E 



A 



£1 \ - / fJa£_ \ B / p a jj 



(13) 



Here A, B and C are the Skyrme parameters which 
are defined according to the ground- state properties of a 
nucleus. Different values of C lead to different equations 
of state. A larger value of C (= 380 MeV) is often dubbed 
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a stiff equation of state. The finite range Yukawa (VV-ufc) 
potential and effective Coulomb potential (Vcoui) read as 

exp{-\n -rj|}//i 



Vy u k 



5>- 

Vcoul — 



ki - r jl/M 



E 

3^3 
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(14) 



(15) 



The Yukawa interaction (with t$= -6.66 MeV and 
p = 1.5 fm) is essential for the surface effects. The 
momentum-dependent interactions (MDI) are obtained 
by parameterizing the momentum dependence of the real 
part of the optical potential. The final form of the po- 
tential reads as follows [12]: 

U MDI « U ln 2 MPi - P 2 ) 2 + l]*(ri - r 2 ). (16) 

where i 4 =1.57 MeV and t 5 =5xl0~ 4 MeV -2 . A param- 
eterized form of the local plus MDI potential is given by 

U = a(£) + + 5\n 2 [e(p/ Po ) 2 ' 3 + l]p/ Po . (17) 

Po Po 

The parameters a, (3, 7, S and e are listed in Ref. [i"2j . 
The momentum-dependent part of the interaction acts 
stro ngly in the cases where the system is mildly excited 
[l5l lilf . In this case, the MDI is reported to generate 
a lot more fragments compared to the static equation of 
state. For a detailed discussion of the different equations 
of state and MDI, the reader is referred to Refs. [3, EH, 
OjJ . The relativistic effect does not play a role in the low 
incident energy of present interest. 

The phase space of the nucleons is stored at several 
time steps. The QMD model does not give any informa- 
tion about the fragments observed at the final stage of 
the reaction. To construct fragments from the present 
phase space, one needs the clusterization algorithms. We 
concentrate here on the minimum spanning tree (MST) 
method and the minimum spanning tree method with 
binding energy check (MSTB) only. 



B. Different clusterization methods 

1. Minimum spanning tree method 

The widely used clusterization algorithm is the MST 
method [13] ■ In the MST method, two nucleons are al- 
lowed to share the same fragment if their centroids are 
closer than a distance r m i n , 



ki-rjl < r„ 



(18) 



where rj and rj are the spatial positions of both nucle- 
ons. The value of r m i n can vary between 2 and 4 fm. 
This method cannot address the question of time scale. 
This method gives a big fragment at high density, which 
splits into several light and medium mass fragments after 
several hundred fm/c. This procedure gives same frag- 
ment pattern for times later than 200 fm/c, but cannot 
be used for earlier times. 
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FIG. 1: (Color online) The time evolution of IMFs (5<A 
<44) for the reaction of 86 Kr+ 93 Nb at incident energy of 75 
MeV/nucleon for different model ingredients . 



Minimum spanning tree method with binding energy 
check 



This is an improved version of the normal MST 
method. First, the simulated phase space is analyzed 
with the MST method and pre-clusters are sorted out. 
Each of the pre-clusters is then subjected to a binding 
energy check [ll], [TtJ : 



N' 



c . in 
Nf 



^ N f 



Vu 



/)] < E, 



bind • 



1=1 j=£i 

(19) 

We take E bind = -4.0 MeV if N f > 3 and E bmd = 0.0 
otherwise. Here N* is the number of nucleons in a frag- 
ment and PStf is the center-of-mass momentum of the 
fragment. This is known as the minimum s pan ning tree 
method with Binding energy check (MSTB) |lil[l7T|. The 
fragments formed with the MSTB are reliable and stable 
at early stages of the reactions. 



III. RESULTS AND DISCUSSION 
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We simulated the central reactions of 
(Ej a6 = 10-55 MeV/nucleon), 40 Ar+ 45 Sc 
35-115 MeV/nucleon), 58 Ni+ 58 Ni (E /ab = 35-105 
MeV/nucleon), 86 Kr+ 93 Nb (E m = 35-95 MeV/nucleon), 
129 Xe+ 118 Sn {E lab = 45-140 MeV/nucleon), 86 Kr+ 197 Au 
(E lab = 35-400 MeV/nucleon) and 197 Au+ 197 Au (E ; „ fc = 
70-130 MeV/nucleon). The energies are guided by exper- 
iments 0-0] ■ For the present study, we use hard (labeled 
Hard), soft (Soft), hard with MDI (HMD), and soft with 
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E c.m (MeV/nucleon) System mass (A) 



FIG. 2: (Color online) The mean IMF multiplicity, (Njmf), 
as a function of incident energy in center-of-mass frame, E c . m . , 
for the reaction 86 Kr+ 93 Nb. Solid (dashed) curves show the 
quadratic fits to the model calculations for MSTB (MST) to 
estimate the peak center-of-mass energy at which the maximal 
IMF emission occurs. 



FIG. 3: (Color online) The E™ 1 (left) and (NiM F } max 
(right) as a function of composite mass of the system (A). 
The different lines in left (right) panels represent the linear fits 
(power-law fits). Comparison of model calculations is made 
with experimental data 0-0] (solid stars). The percolation 
calculations [(J (solid triangles) are also shown in figure. 



MDI (SMD) equations of state. We also use the stan- 
dard energy-dependent Cugnon cross section (cr( l 7 ^ e ) [l6[ 
and constant isotropic cross section of 55 mb strength in 
addition to two different widths of Gaussian L = 1.08 
and 2.16 fm 2 (L broad ). The superscripts represent cross 
section. The phase space is clusterized using the cluster- 
ization methods described previously. The reactions are 
followed until 200 fm/c but the conclusions do not change 
when the reaction is complete, employing the validity of 
both algorithms. 

In fig. 1, we display the time evolution of IMFs 
for the reaction 86 Kr+ 93 Nb at incident energy of 75 
MeV/nucleon employing the MST method. In fig. 1(a), 
we display the model calculations using Hard c " s (solid 
line) and Soft Cua (dashed line). From fig. 1(a), we see 
that the number of IMFs is larger in case of Soft com- 
pared to that of Hard. This is because of the fact that soft 
matter is easily compressed, resulting in greater achieved 
density, which in turn leads to the large number of IMFs 
compared to that in the Hard. It is worth mentioning 
here that the effect could be opposite at higher ener- 
gies, because at higher energies the IMFs may further 
break into LCPs and free nucleons. In fig. 1(b), we 
display the results for Hard c " 9 and Hard 55 (dotted line) . 
As evident from the fig. 1(b), the number of IMFs is 
nearly same for both types of cross sections. This may 



be because, for the central collisions, since the excita- 
tion energy is already high, different cross sections have 
a negligible role to play. In fig. 1(c), we display the re- 
sults for the Hard case along with two different widths of 
the Gaussian,that is, L and \} road (dash-dotted line). We 
find that the width of Gaussian has a considerable impact 
on fragmentation. As we change the Gaussian width (L) 
from 4.33 to 8.66 fm 2 , the multiplicity of IMFs is reduced 
by » 30%. Interestin gly, the kaon yield also gets reduced 
by the same amount [18j |. Owing to its large interaction 
range, an extended wave packet (i.e. L fcroad ) connects a 
large number of nucleons in a fragment, as a result, it gen- 
erates heavier fragments compared to what is obtained 
with a smaller width. It is worth mentioning here that 
the width of the Gaussian has a considerable effect on 
the collective flow [HI, [l9| as well as on pion production 
[H|,[2(|. In fig. 1(d), we display the results using Hard 
and HMD (dash-dot-dotted line). Again the number of 
IMFs are nearly same for both equations of state (EOS). 
This is expected because the effect of MDI is small at 
these energies. However, the scenario is completely dif- 
ferent at high energies; at high energies, owing to the 
repulsive nature of MDI, there is a large destruction of 
initial correlations and the additional momentum depen- 
dence further destroys the correlations reducing further 



5 



the multiplicity of IMFs. This leads to the emission of 
lots of nucleons and LCPs [2l| . 

In fig. 2, we display the average multiplicity of IMFs, 
(N/mf), as a function of incident energy in the center- 
of-mass frame (E c m ) for the 86 Kr+ 93 Nb reaction em- 
ploying the MST (open symbols) and the MSTB (solid 
symbols) methods. Figures. 2(a) and 2(b) are for Soft 9 
and Hard " 9 , respectively. Lines represent the quadratic 
fit to the model calculations. In both cases, the num- 
ber of IMFs first increases with incident energy, attains 
a maximum, and then decreases, in agreement with pre- 
vious studies 0-0, M, ED]- Clearly, (N/ MF ) is greater for 
the MST method than for MSTB, because in the case 
of the MSTB, along with spatial correlations, an addi- 
tional check for binding energy is also used; therefore, 
it filters out the loosely bound fragments which decay 
later. Hence, the fragments obtained with THE MSTB 
are properly bound. A similar trend is obtained for all 
other reactions as well as different model ingredients used 
in the present study but is less pronounced in lighter 
systems like 20 Ne+ 20 Ne and 40 Ar+ 45 Sc as compared to 
heavier systems. However, for the Gaussian width L broad , 
the value of (N/ mf ) is nearly zero in this incident energy 
range using THE MSTB (not shown here). This is be- 
cause an extended wave packet (i.e. L broa<i ) connects a 
large number of nucleons in a fragment; as a result, it 
generates heavier fragments and the additional binding 
energy check further excludes the unbound fragments. 

In fig. 3, we display the peak center-of-mass en- 
ergy E™^f (left panels) and peak multiplicity of IMFs 
(N/mf)™" 1 (right panels) as a function of the combined 
mass of the system employing the MST method. In the 
left panels, lines represent linear fitting proportional to 
mA and in the right panels, lines represent power-law fit- 
ting proportional to A r . The multiplicity of IMFs, in the 
case of 20 Ne+ 20 Ne and 40 Ar+ 45 Sc, is obtained by exclud- 
ing the largest and second largest fragment, respectively, 
to get the accurate information about the system size de- 
pendence. (NiMF) max and corresponding E™^ are ob- 
tained by making a quadratic fit to the model calculations 
for (N/mf) as a function of (E c m ). From the left panels, 
we find that the mass dependence of E" l ^f is insensitive 
to different EOS (fig. 3a), nn cross section (fig. 3b), as 
well as the width of the Gaussian [fig. 3c]. It is slightly 
sensitive to MDI because, for heavy systems, the value of 
^c.m. i s greater, as a result of which is that the effect of 
MDI becomes non- negligible and creates the (^iMF) max 
at smaller energies. From the right panels where we dis- 
play the mass dependence of (NiMF) max , we find that 
the peak multiplicity is insensitive to cross section (fig. 
3d) and MDI (fig. 3h) (for explanation see discussion of 
fig. 1). It is slightly sensitive to the EOS (fig. 3b) but 
highly sensitive to the Gaussian width (fig. 3f). On in- 
creasing the width of the Gaussian, (Njmf)™ 021 reduces 
to a large extent. As discussed earlier, an extended wave 
packet (i.e. L fcroad ) will connect a large number of nucle- 
ons in a fragment; as a result, it generates heavier frag- 
ments compared to what one obtains with smaller width. 
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FIG. 4: (Color online) The E££* (left panels) and {N IM F) max 
(right panels) as a function of composite mass of the system 
(A) using Soft " 9 (upper panels) and Hard c " 9 (lower panels) 
employing THE MSTB and THE MST methods. Lines have 
same meaning as in fig. 3. Comparison of model calculations 
is made with experimental data 043] (solid stars). 



From fig. 3, we see that E™^f shows linear dependence 
(oc mA) whereas (N;mf)' mi follows power-law behavior 
(oc A r ) with t nearly equal to unity. 

In fig. 3, the model calculations are also compared with 
experimental data 0-0] • It is clear from fig. 3 that model 
calculations for E™^ agree with experimental data 0-0] 
whereas in the case of {NiMF) max , as the system mass 
increases, the difference between model calculations and 
experimental results continues to increase. This behavior 
is consistent for all the different choices of model ingredi- 
ents. This may be because the fragments obtained with 
the MST method are not reliable because this method 
makes sense only when matter is diluted and well sep- 
arated. This is true only in the case of high beam en- 
ergy and in central collisions. Therefore, we have to look 
for other methods of clusterization. As mentioned ear- 
lier, the fragments obtained with the MSTB are properly 
bound and reliable. So, as a next step, we check the sys- 
tem size dependence of E™^ x and (Njjhf) 1 ™" 1 by using 
the MSTB for clusterization. 

In fig. 4, we display the E™^ x (left panels) and 
(N/MF) max (right panels) for Soft c " s (upper panels) and 
Hard Ctl9 (bottom panels) as a function of the combined 
mass of the system. Solid (open) symbols represent the 
MSTB (MST method). From the left-hand panels, we 
find that E™£f remains insensitive to the choice of clus- 
terization method. The same is true for (N]MF) max 
(right-hand panels) but in low-mass region. As the sys- 
tem mass increases, (N/mf)™" 1 becomes more and more 
sensitive to the method of clusterizaton. The MSTB ex- 
cludes the loosely bound fragments, thus reducing the 
peak IMF multiplicity. The effect is uniform both for the 
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FIG. 5: (Color online) The largest fragment and multiplic- 
ities of free-nucleons, LCPs, MMFs, HMFs, and IMFs as a 
function of composite mass of the colliding nuclei (A) for dif- 
ferent reactions at their respective E™^ 1 (solid circles). Lines 
represent the power-law fits (oc A T ). 



EOS and for different cross sections (not shown here). 

In fig. 5, we display peak multiplicity (obtained by em- 
ploying the MSTB) as a function of composite mass of 
the system for various fragments consisting of the largest 
fragment (A max ) (fig. 5a), free-nucleons (1<A <1) (fig. 
5b), light charged particles (LCPs) (2<A <4) (fig. 5c), 
medium mass fragments (MMFs) (5<A <9) (fig. 5d), 
heavy mass fragments (HMFs) (10<A <44) (fig. 5e) and 
intermediate mass fragments (IMFs) (5<A <44) (fig. 5f) 
for Hard " 9 . Lines represent the power-law fitting pro- 
portional to A T . Interestingly, the peak multiplicities 
of different fragments follow a power-law (oc A T ). The 
power-law factor r is almost unity in all cases except 
pmax j? or wri ich th ere jg n0 clear system size dependence. 
The system size dependence of various fragments was also 
predicted by Dhawan and Puri [22|]. Their calculations 
at the energy of vanishing flow (i.e., the energy at which 
the transverse flow vanishes) clearly suggested the exis- 
tence of a power-law system mass dependence for various 
fragment multiplicities. 

To check the stability of fragments, we display in 
fig. 6, the binding energy per nucleon as a function of 




Time (fm/c) 



FIG. 6: (Color online) The time evolution of binding en- 
ergy per nucleon for LCPs, MMFs, HMFs, and IMFs. Here 
reactions of 20 Ne+ 20 Ne, 40 Ar+ 45 Sc, 58 Ni+ 5S Ni, 86 Kr+ 93 Nb, 
129 Xe+ 118 Sn, 86 Kr+ 197 Au, and 197 Au+ 197 Au are simulated 
for central geometry at their corresponding E™^ x . 



time for LCPs, MMFs, HMFs, and IMFs. The reac- 
tions of 20 Ne+ 20 Ne, 40 Ar+ 45 Sc, 58 Ni+ 58 Ni, 86 Kr+ 93 Nb, 

129 Xe+ 118 Sn; 86 Kr+ 197 Au and 197 Au+ 197 Au &re simu _ 

lated at laboratory energies corresponding to their E™^ 1 
values, which are approximately 24, 46, 69, 78, 96, 124, 
and 105 MeV/nucleon, respectively. We find that, even 
at 200 fm/c, small fragments are still not cold and take a 
very long time to cool down, whereas the heavy fragments 
are properly bound, having binding energy per nucleon 
around -5 to -7 MeV. 



IV. SUMMARY 

We simulated the central reactions of nearly symmet- 
ric, and asymmetric systems over the entire periodic table 
at different incident energies for different EOS, nn cross 
sections, and different widths of Gaussians. We observed 
that the multiplicity of IMFs (3 <Z< 20) shows a rise and 
fall with increase in beam energy in the center-of-mass 
frame, as already predicted experimentally and theoret- 
ically. We also studied the system size dependence of 
peak center-of-mass energy EJT^f and peak IMF multi- 
plicity (N/mf)" 1 " 1 - It was observed that E™^f increases 
linearly with system mass whereas a power-law (oc A r ) 
dependence was observed for (N/mf)" 1 " 1 with r ~1.0. 
We compared the system size dependence of E™^ and 
(NiMF) max for MST and MSTB methods and found that 
MSTB reduces the (N/mf)™" 1 especially in heavy sys- 
tems because loosely bound fragments get excluded in 
MSTB. The power-law dependence is also observed for 
fragments of different sizes at the energy for which the 
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production of IMFs is a maximum and the power-law 
parameter r is found to be close to unity in all cases ex- 
cept A max . The stability of fragments is also checked 
through binding energy per nucleon. We observed that, 
at 200 fm/c, small fragments are still not cold and they 
take a very long time to cool down, whereas the heavy 
fragments are properly bound. 
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